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The effects of astrophysical uncertainties on the exclusion limits at dark matter direct detection 
experiments are investigated for three scenarios: elastic, momentum dependent and inelastically 
scattering dark matter. We find that varying the dark matter galactic escape velocity and the 
Sun's circular velocity can lead to significant variations in the exclusion limits for light (< 10 
GeV) elastic and inelastic scattering dark matter. We also calculate the limits using one hundred 
velocity distributions extracted from the Via Lactea II and GHALO N-body simulations and find 
that a Maxwell-Boltzmann distribution with the same astrophysical parameters generally sets less 
constraining limits. The elastic and momentum dependent limits remain robust for masses > 50 
GeV under variations of the astrophysical parameters and the form of the velocity distribution. 
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I. INTRODUCTION 



Identifying the dark matter in our Universe is one of 
the great challenges facing both theorists and experi- 
mentalists. Dark matter direct detection experiments 
are now probing dark matter-nucleon cross sections ex- 
pected of Weakly Interacting Massive Particles (WIMPs) 
and the results from the CDMS II [1^ and CoGeNT [2] 
experiments may be hinting that a conclusive detection 
of dark matter is within reach of the next generation 
of experiments. The DAMA collaboration [3] have also 
continued to measure an annual modulation signal, ob- 
served over the past thirteen years with an 8.9(t signifi- 
cance which has many of the properties expected from a 
dark matter origin. Although now inconsistent with the 
null results from other experiments assuming a vanilla 
elastically scattering WIMP, this persistent signal has 
inspired many models [1U12) which can give rise to dra- 
matically different experimental signals from the usual 
assumptions of elastically scattering spin-independent or 
spin-dependent interactions. Although these models may 
be capable of reconciling the DAMA result with other ex- 
periments, we consider them to be interesting, indepen- 
dent of the DAMA result, since they serve to highlight 
the different possible signals that could be produced at 
dark matter direct detection experiments. 

As well as depending on particle physics models, the 
recoil spectrum and scattering rate depend on the as- 
trophysical parameters and the form of the local veloc- 
ity distribution the particle dark matter in the galactic 
halo takes. Unfortunately, these parameters still suffer 
from large uncertainties, while high resolution numeri- 
cal N-body simulations have shown that the velocity dis- 
tribution deviates significantly from the default choice, 
namely, a Maxwell-Boltzmann distribution. 

Given this, in this paper we investigate the effect 
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of astrophysical uncertainties on three scenarios: elas- 
tic spin-independent scattering; and, as proxies for the 
many varied models, inelastic and momentum depen- 
dent (aka form factor) spin-independent scattering. We 
vary the galactic escape velocity, the Sun's velocity rel- 
ative to the galactic centre, and the local dark matter 
density within their experimental errors. We then com- 
pare the exclusion limits for these scenarios assuming a 
Maxwell-Boltzmann velocity distribution with those ex- 
tracted from two of the highest resolution dark matter N- 
body simulations. Via Lactea II and GHALO. We partic- 
ularly focus on four experiments with different target nu- 
clei which span a large range of atomic number A to com- 
pare the relative change at different experiments: CDMS 
Il-Si; CDMS II-Ge; XENONIO; and CRESST-II, with 
silicon {A = 28), germanium {A = 72), xenon {A = 131) 
and tungsten [A = 184) target nuclei respectively.^ 

In Fig. [T] we have plotted the differential scattering 
rate at a germanium detector to illustrate the differences 
between the three scenarios. Inelastic scattering of the 
dark matter to a state ~ 100 keV heavier [4] changes 
the phenomenology in three principal ways: it leads to 
a suppression of scattering rates at low mass target nu- 
clei experiments; an increase in the ratio of the modu- 
lated to unmodulated scattering rate; and a kinematical 
suppression of the recoil spectrum at low recoil energies. 
Although often invoked as a way of reconciling the modu- 
lation signal seen by DAMA with the null results at other 
experiments |15H19) , the latest (preliminary) results from 
CRESST-II [20] tentatively rule out this interpretation. 
However, having this splitting is natural in many models 
[2TH30] and leads to much other interesting phenomenol- 
ogy, for example through capture in astrophysical bodies 
[?m35| . Therefore, we emphasise that even if inelastic 



^ Since we are seeking to understand the effects of the astrophysical 
uncertainties, we do not specifically address whether astrophysical 
uncertainties can improve the fit between the CoGeNT and DAMA 
signals. For analyses in this direction, see |13II14| . 
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FIG. 1: (Colour online). The differential scattering rate at a 
germanium detector for M^ — 100 GeV for elastic, inelastic 
and momentum dependent (MD) spin-independent scatter- 
ing. The normahsation of each curve is different. 



scattering can not reconcile the DAMA result with other 
experiments, it remains an interesting and viable possi- 
bility to explore. 

In momentum dependent scattering [6l [7] , the differ- 
ential cross section is suppressed by additional powers of 
g2 = 2MNEji, leading to a suppression of the exclusion 
limits from experiments which rely on having a sensitivity 
at low recoil energies, in our case XENONIO and CDMS 
Il-Si. This dependence can arise for example through 
dipole or charge radius interactions [36 . 

In Section |ll] we briefly review the theory required to 
calculate scattering rates, then in Section [Till we review 
our current knowledge of the relevant astrophysical pa- 
rameters and demonstrate and explain their effect on the 
experimental exclusion curves for elastic, inelastic and 
momentum dependent scattering. In Section |IV| we ex- 
amine the effect of using velocity distributions from the 
Via Lactea II and GHALO simulations. Finally, Ap- 
pendix |A] summarises the relevant experimental details, 
where in particular, we call attention to the reanalysed 
early CDMS II runs and discuss uncertainties affecting 
the XENONIO exclusion limits at low masses. 



II. REVIEW OF STANDARD FORMALISM 

In this section we briefly review the theory needed to 
calculate the number of events expected at a given ex- 
periment. The differential scattering rate is given by 
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where Nt is the number of target nuclei, p^ is the lo- 
cal dark matter density, My^ is the dark matter mass, 
V = |i/|, and f^{v,Ve) is the local dark matter velocity 
distribution in the detector rest frame. For a nucleus 
to recoil with an energy En, the incident dark matter 



particle must have a minimum speed given by: 
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Here Mjv is the nucleus mass, p,^ the dark matter- 
nucleus reduced mass, and we include (5, the mass dif- 
ference between the incoming and outgoing dark matter 
particle, to allow for inelastic scattering. 

We will concentrate on coherent spin-independent in- 
teractions in the remainder of this paper, in which case, 
the dark matter-nucleus differential cross section is given 
by 
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where cr„ is the dark matter-neutron cross section at zero 
momentum transfer in the elastic limit, fine is the dark 
matter-nucleon reduced mass and /„ and fp parameterize 
the relative dark matter coupling to neutrons and pro- 
tons respectively. For simplicity, we shall assume /„ = fp 
for the remainder of this paper, but one should bear in 
mind that in specific particle physics models, this may 
not be the case. Fortunately the rescaling of the cross 
section for different values of /„ and fp is trivial. F{Er) 
is the nucleus form factor which for tungsten and xenon 
nuclei, we take as the Fermi Two-Parameter distribution 
with parameters taken from the Appendix of ^37j . For 
germanium and silicon we use the Lewin-Smith [38] pa- 
rameterisation of the Helm form factor: 
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where q = ^/2M]sEr, R = J c^ + ^-k'^ a^ - bs'^ , c = 

1.23^1/3 - 0.60 fm, s = 0.9 fm and a = 0.52 fm. We 
find that this parameterisation agrees well (to within 2% 
or better) with the 'spherical Bessel expansion' and 'sum 
of Gaussian expansion' from |37j , in stark contrast to the 
commonly used values R = y/r"^ — 5s^, s = 1 fm and 
r = l.2A^/^ which overestimates F'^(Er) by up to 20% 
at Er = 100 keV. As pointed out in [3T, the uncertainty 
in the exclusion limits due to the form factor parameter- 
isation can be significantly reduced by a judicious choice, 
leaving, as we will show, uncertainties in the astrophysi- 
cal parameters as the dominant source of uncertainty in 
the exclusion limits. 

Finally, to allow for momentum dependent scattering, 
we follow j^ by parameterising the momentum depen- 
dence by replacing the differential cross section with 



(5) 



where g^e/ — 100 MeV. Here we only consider the case 
n — 1, although models with n — 2 have also been con- 
structed. 




A. Astrophysical formalism 

The simplest model of the dark matter distribution 
is the so called 'standard halo model' (SHM), which as- 
sumes an isotropic, isothermal sphere for the dark matter 
distribution, and leads to a Maxwell-Boltzmann velocity 
distribution in the galactic frame, truncated at the galaxy 
escape velocity Vesc'- 
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Here A'' is a normalisation constant and vq is the circular 
speed of the Sun about the centre of the galaxy (which 
throughout this paper we assume is the same as the local 
standard of rest). High resolution N-body simulations of 
the galactic dark matter structure have shown that the 
dark matter velocity distribution departs from the SHM 
[391 - I4T] . and the effects can be large, especially in scenar- 
ios such as inelastic dark matter where only the high ve- 
locity tail of the velocity distribution is sampled [T71 WH - 
H5] . Therefore in Section IV we explore how the exclusion 
limits vary when we use data from the Via Lactea H [50] 
and GHALO [51 simulations, two of the highest resolu- 
tion simulations of galactic dark matter structure. 

The velocity distribution entering Eqn. [l] is evalu- 
ated in the frame of the detector, therefore we need to 
transform the galactic velocity distibution (Eqn. [6]) to 
the Earth frame. f^{v,Ve) and f{v) are simply related 
through: 
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where Ve = vq+v^- vq — wo + %) is the sum of the Sun's 
circular velocity vq relative to the galactic centre, and 
peculiar velocity %> relative to Sun's circular velocity. 
v^ is the Earth's velocity relative to the Sun's rest frame 
which we calculate following the prescription in Appendix 
B of Ref. [38. Combining Eqns. Ill Island [tI the relevant 
quantity which enters into the integration in dR/dEn is: 
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In Appendix [B] we present an analytic expression for 
C{Er) for the SHM velocity distribution, as defined by 
Eqn. [6] which to the best of our knowledge has not pre- 
viously been presented in the literature. 



III. ASTROPHYSICAL UNCERTAINTIES 

The quantities entering into this calculation from as- 
trophysics are: p^, the local dark matter density, with 
canonical value p^ = 0.3 GeV/cm'^; the Sun's circular 
velocity with respect to the galactic rest frame with fidu- 
cial value Vq = (0, 220, 0) km/s; the Sun's peculiar veloc- 
ity relative to the Sun's circular velocity, which is usually 



given by u® = (10.0, 5.25, 7.17) km/s (55]; and the escape 
velocity Vesc-, which is often chosen in the range from 500 
km/s to 650 km/s. In this section we will discuss recent 
determinations of these parameters and investigate the 
changes in the exclusion limits when we vary them within 
their allowed ranges, for elastic, momentum dependent 
and inelastic scattering dark matter with 5 = 130 keV. 



A. The galactic escape velocity 

The RAVE survey |S3] has determined a 90% con- 
fidence interval for the galactic escape velocity of 
498 km/s < Vesc < 608 km/s with a median likelihood 
of 544 km/s. In Fig. l2] we plot the exclusion limits for 
CDMS Il-Si in brown, CDMS II-Ge in blue, CRESST- 
II in green and XENONIO in red for elastically scatter- 
ing (left and right upper panels), momentum dependent 
(lower left panel) and inelastically scattering dark mat- 
ter with a mass splitting 5 = 130 keV (lower right panel) 
for three values of Vesc', 498 km/s (dotted); 544 km/s 
(solid); and 608 km/s (dashed) while keeping p^ — 0.3 
GeV/cm'^ and Vq = 220 km/s fixed. For clarity, we have 
only shown the limits for CDMS Il-Si and XENONIO 
at low masses (upper left panel), since the CDMS II-Ge 
limits behave similarly to those from XENONIO in this 
mass range. CDMS Il-Si is not able to set any limits 
for inelastic dark matter with 6 = 130 keV since its low 
mass means the minimum speed required to scatter off a 
nucleus is higher than the galactic escape velocity. 

The limits for momentum dependent and elastically 
scattering dark matter are very insensitive to Vesc for all 
experiments and it is only for masses M^ < 10 GeV that 
the three lines can be separately resolved. In compar- 
ison, inelastically scattering dark matter is much more 
sensitive, with two principal features to note: the curves 
shift left (right) and down (up) when Vesc is increased 
(decreased). 

In Fig. [3] the fractional change in an relative to the 
value at Vesc = 544 km/s is plotted as a function of 
mass for the inelastically scattering case. The dotted 
and dashed lines show the ratio for Vgsc = 498 km/s and 
608 km/s respectively. We see that the XENONIO and 
CRESST-H limits change in a similar way, varying by 
< 20% at high masses, and remaining fairly constant 
until ~ 100 GeV, when a very rapid increase occurs. 
In comparison, decreasing the escape speed dramatically 
weakens the CDMS II-Ge limits, by 80% at high masses 
(M^ ~ 1000 GeV) and larger at lower masses. Increasing 
Vesc has a less dramatic effect with a change of ^ 40% 
for CDMS H-Ge until ^ 150 GeV when a rapid increase 
happens. These effects are simple to understand qualita- 
tively and quantitively. 

We will first explain why the curves move vertically in 
Fig. [2j If there is a lower escape velocity, there will be 
fewer particles in the halo capable of scattering, there- 
fore a larger cr„ is required to produce the same number 
of events at a given experiment. Since the Maxwell- 
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FIG. 2: (Colour online). Varying Vesc within the 90% confidence limit found by the RAVE survey: Vesc = 498 km/s (dotted), 
544 km/s (solid) and 608 km/s (dashed). Varying Vesc has a negligible effect on the hmits for momentum dependent and 
elastically scattering dark matter except at low masses, where the limits are shifted horizontally by ~ 0.5 GeV. For inelastically 
scattering dark matter, the effect is much larger, especially for CDMS II-Ge. The two principal features to note are a shift left 
(right) and down (up) when Vesc is increased (decreased). 



Boltzmann distribution is exponentially suppressed at 
high velocities, decreasing Vesc will have a very small ef- 
fect in the number of particles able to scatter, unless Vmin 
is close to "yesci in which case the fractional change can be 
substantial. For elastic and momentum dependent scat- 
tering dark matter, this occurs at M^ < 10 GeV, and is 
true for all masses when 5 = 130 keV. CDMS II-Ge is par- 
ticularly affected relative to XENONIO and CRESST-II 
because its germanium target is much lighter than xenon 
and tungsten, so Vmm is higher, and much closer to Vesc- 
Similarly if v^sc is raised, there are extra particles in the 



halo which can scatter, so a lower (t„ will produce the 
same number of events in an experiment. Due to the 
exponential suppression of the Maxwell-Boltzmann dis- 
tribution, increasing v^g^ has a smaller affect than de- 
creasing it by a similar amount. 

Unfortunately, the intricacies of the 'pmax' method 
|54) for setting exclusion limits prevents very accurate 
estimates. However, we can give a quantitive estimate 
which should serve as a useful guide to the size of shift 
expected when v^.^^ is varied. Since the combination (7„C 
enters dR/dEji (from Eqns.lslandls]), the vertical shift is 
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FIG. 3: (Colour online). The fractional change in (t„ rela- 
tive to value at Veac = 544 km/s for Vesc = 498 km/s (dotted) 
and v^sc ~ 608 km/s (dashed) for inelastically scattering dark 
matter with 5 = 130 keV. The blue, red and green lines show 
the change for CDMS II-Ge, XENONIO and CRESST-II re- 
spectively. 



given by 
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where an analytic expression for ^ is given by Eqn. |B5[ 
For example, for CRESST-II with v^^n from Eqn. [2| 
Mx = 500 GeV, S = 130 keV, and Er = 25 keV? 
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good agreement with the values 0.15 and —0.08 pre- 
sented in Fig. |3] In comparison for CDMS II-Ge, with 
M^ = 500 GeV, 6 = 130 keV and Er = 80 keV 2, we 



nnd Act.,,, , — 0.75 and Act 



= -0.29. This 
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compares to the values of 1.1 and —0.39 shown in Fig. [3j 
a difference of around 40%. 

To understand the horizontal shift, we will consider the 
energy of the recoiling nucleus, which is ultimately what 
experiments measure. In the Earth frame, the nuclear 
recoil energy is given by 
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where v^ is the dark matter speed in the Earth's frame 
and Or is the nuclear recoil angle. If the recoil energy is 
the same, we expect a similar signal at an experiment. 
Increasing Vesc will increase the average speed of dark 
matter particles in the halo, but this increase will only 
be perceptible if Wmin is close to Vgsc, as for inelastic and 
light dark matter. When the average speed increases, the 
average kinetic energy will be the same for a lighter dark 
matter particle, and since we expect the nuclear recoil 
energy to depend on a simple way on the dark matter's 



initial kinetic energy, the exclusion curves shift to the 
left. 3 

Quantitively, the mass shift is found from differentiat- 
ing Eqn. [lOl leading to 
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where AM^ = M^2 — M-^i, and Av^ = v^2 — Vx^ '^o'"" 
responds to the shift in the average speed of the dark 
matter responsible for measurable events. 

For CDMS H-Ge with 6 = 130 keV, M^ = 90 GeV, 
the only particles which can scatter have a speed v^ ~ 
Vesc + Ve — 800 km/s. Changing v^sc from 544 km/s by 
±64 km/s leads to Au^ = ±64 km/s, which from Eqn. 
[TT| gives AM^ ~ =f17 GeV, in reasonable agreement with 
what we find in Fig. [2J For the elastic case Aw-^ '-^ 
because the average speed doesn't depend sensitively on 
the escape speed, unless Vmm ~ Vesc + Ve which occurs 
at M^ - 10 GeV. For XENONIO at M^ = 7 GeV, 5 = 
keV, v^ = 800 km/s and Av^ = ±64 km/s, this leads to 
AM^ ~ =F0.6 GeV, again in reasonable agreement with 
Fig. [I 



B. The Sun's peculiar velocity 

The pecuhar velocity w® = {U, V, W^)® is often as- 
sumed to be well known: i/@ = (10.0,5.25,7.17) km/s 
with small errors ^ 0.5 km/s [55]. However three recent 
papers have suggested that V® may have been underesti- 
mated by -- 7 km/s [57H5^ . Ref.[^ found that the clas- 
sical determination in |52j underestimated V® by ignoring 
the metallicity gradient in the Milky Way disc. Further- 
more when w® is increased, [S7| and '55] found better fits 
to the observations of masers in high mass star forming 
regions and the Milky Way distribution functions respec- 
tively. Therefore unless otherwise stated, in the remain- 
der of this paper we will use this newly determined value 
for the Sun's peculiar velocity: w® = (11.1,12.24,7.25) 
km/s with uncertainties ~ (1,2,0.5) km/s. 



C. The Sun's circular velocity 

The Sun's circular velocity can be determined in a 
number of ways, which give slightly different values from 
the fiducial value of vq = 220 km/s. The ratio of the 
total velocity of the Sun about the galactic centre to the 
distance of the Sun from the galactic centre has been de- 
termined from the apparent proper motion of Sgr A* to 
hevQ/Ro = {vq + V(^)/Ro = 30.2±0.2 km s^^kpc^i [50]. 
Unfortunately Rq is still poorly determined with the most 
recent studies of stellar orbits obtaining Rq = 8.4 ± 0.4 



^The peak energy of the recoil spectrum dR/dEji |16| 



^A similar argument is presented in 1551 1561 
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FIG. 4: (Colour online). The exclusion limits for va = 195 km/s (dotted), iiq = 220 km/s (solid) and wo = 255 km/s (dashed). 
Varying wo leads to a larger vertical shift in the exclusion limits compared to varying v^sc for all panels (cf. Fig. [2|. The 
horizontal shift for all experiments at masses ~ 10 GeV for elastic and momentum dependent scattering is ~ 1 — 2 GeV, larger 
than when Vesc is varied (cf. Fig.[2|. In contrast, when varying vq in the inelastic case, the horizontal shift is smaller compared 
to varying v^ac- 



kpc f5P and i?o = 8.33 ± 0.35 kpc f5?. Combining these 
results, and using the new value of V® = 12.24 km/s, we 
find I'o = 242 ± 12 km/s and wq — 239 ±11 km/s respec- 
tively. This is in comparison to the analysis of the GD-1 
stellar stream which found vq = 221 ± 18 km/s [63_, and 
the best fit to masers in the high mass forming regions 
which found the range t^o = 225 ± 29 km/s [STj. Wary 
of the possibility of unknown systematic errors affecting 
one of these measurements of uq, we take a conservative 
approach by giving each the same weight and hence in 
Fig. [4] we investigate how the exclusion limits change for 



three values of vq\ 195 km/s (dotted); 220 km/s (solid); 
and 255 km/s (dashed), while keeping p^ = 0.3 GeV/cm^ 
and Wesc = 544 km/s fixed. Again the left and right up- 
per panels shows elastically scattering dark matter while 
the lower left and right panels show momentum depen- 
dent and inelastically scattering dark matter for b = 130 
keV. We see that varying wg has more of an effect on a^ 
than changing v^sc for all cases (cf. Fig. [2]), and that once 
again there are two principal features to note: increasing 
(decreasing) wq causes the exclusion curves to shift down 
(up) and left (right). 
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FIG. 5: (Colour online). The fractional change in (j„ rel- 
ative to value at vo — 220 km/s for vq = 195 km/s (dot- 
ted) and Wo = 255 km/s (dashed) for elastically (top panel), 
momentum dependent (middle panel) and inelastically scat- 
tering dark matter with S — 130 keV (bottom panel). The 
blue, brown, red and green lines show the change for CDMS 
II-Ge, CDMS Il-Si, XENONIO and CRESST-II respectively. 
For clarity we only show the change for CDMS II-Ge and 
XENONIO in the top panel and CDMS II-Ge and CDMS II- 
Si in the middle panel. 



Changing vq has two effects: vq is needed to boost 
from the galactic frame to the Earth frame (Eqn. [7]) , so 
increasing vq means there are more particles at speeds 
above Vmin available to scatter in the Earth's frame; and 
in the SHM, vq determines the dark matter dispersion 

velocity (cTdis = a/I'^o): so increasing vq also increases 
the range of the speed of particles in the halo, with more 



Vesc, as for inelastic dark matter and light dark matter, 
changing vq has a similar effect to changing Vesc since it 
also increases or decreases the number of particles which 
can scatter. However since we are not just changing the 
tail of an exponential distribution, but rather its overall 
shape and position, the effect of changing vq is larger than 
varying v^sc- This is demonstrated in the bottom panel 
of Fig. [5] where we have plotted the fractional change in 
cr„ relative to the value at vq ~ 220 km/s for vq = 195 
km/s (dotted) and vq = 255 km/s (dashed) assuming 
6 — 130 keV. Comparing this to Fig. [3] we see that it has 
the same features, but the effects are exaggerated. 

The top and middle panel of Fig.|5]shows the fractional 
change in an relative to the value at vq = 220 km/s for 
Vq = 195 km/s (dotted) and vq — 255 km/s (dashed), 
assuming elastic and momentum dependent scattering. 
For clarity, in the top panel, we have only plotted the 
fractional change for CDMS II-Ge (blue) and XENONIO 
(red) since the fractional changes for CRESST-II and 
CDMS Il-Si are similar to those of XENONIO. In the 
middle panel, we have only plotted the change for CDMS 
II-Ge and CDMS Il-Si. In the top panel, for masses > 75 
GeV, the fractional change in i7„ when increasing or de- 
creasing Wo remains at ~ 10% for all experiments. How- 
ever, for lower masses, the fractional change at all ex- 
periments increases dramatically, reaching ~ 100% at 10 
GeV. This sudden increase is caused by the horizontal 
shift in the exclusion limits which occurs for the same 
reasons as outlined in Section IIII Al above. In the mid- 
dle panel, the CDMS II-Ge curves are similar to those 
in the top panel (although not shown, XENONIO and 
CRESST-II are also similar), but those for CDMS Il-Si 
are quite different. This is because for momentum depen- 
dent scattering, the events at higher recoil energies are 
more important. To scatter with a larger recoil energy, a 
larger speed is required. For CDMS Il-Si, this minimum 
speed is close to Vesc, so the limits are effected more when 
Vq is varied, compared with elastic scattering. 

We can use the formulas from Section PlI Al to estimate 
the vertical and horizontal shifts of the exclusion curves 
for inelastic scattering. In particular for XENONIO at 
S = 130 keV and Er = 35 keV,l2l we 

0.56, in good 

agreement with the values 1.32 and —0.60 presented in 
the bottom panel of Fig. [5j while for the horizontal shift 
at M^ = 65 GeV, S = 130 keV, Av^ = ±35 km/s and 
v^ ~ 700 km/s we find AM^ ~ t5 GeV. For XENONIO 
at M^ = 12 GcV, 6 = keV, Av^ = ±35 km/s and 
v^ ~ 500 km/s we find AM^ ~ =f1 GeV. These estimates 
are all in reasonable agreement with what is observed in 

Fig.m 



Mx = 400 GeV 

n 1 A 195km/s -, rn j A 2bdkm/s 

find Act22o^,„^^ = 1-59 and Aa^^^^J^^ = 



at higher and lower speeds. Therefore if Vmin is close to 



D. The local dark matter density 

We next consider the local dark matter density p^^. 
Unfortunately, this is the least well known astrophysical 
parameter with an often quoted uncertainty of a factor 



of 2 or 3 in the fiducial value p^ = 0.3 GeV/cm^ [Ml 155] , 
Two recent studies using different techniques found p^ = 
0.3 ± 0.1 GeV/cm3 [55] and p^ = 0.43 ± 0.15 GeV/cm^ 
[57] . consistent with p^ — 0.3 GeV/cm'^ with an error 
of a factor of 2. Dark matter N-body simulations from 
the Aquarius Project have shown that p^^ should be very 
smooth at the Sun's position, varying by less than 15% at 
the 99.9% confidence level from the average value over an 
ellipsoidal shell at the Sun's position [iO]. Therefore we 
can be reasonably confident that the Earth is not sitting 
in a particularly over or under dense region of the halo. 
Since the combination p-^Un enters into dR/dEj^ (Eqn. 
fl]), the error in p^ corresponds directly to an uncertainty 
of a factor of 2 in the dark matter-neutron cross section 
cr„, affecting all experiments equally. 



E. Discussion: The correlation between vo and p^ 

We have followed the standard procedure of varying 
these astrophysical parameters independently, however, 
Refs. [57| [5D] found that there is a strong correlation 
between the Sun's circular speed vq and the distance from 
the galactic centre Rq: vq/Ro w 29 kms~^kpc~^. Since 
p^ also depends on Rq, a more careful approach should 
take into account the change in p^ as we vary vq. Here 
we estimate this correction. 

Assuming vq = 230 km/s and p^ — 0.3 GeV/cm'^ at 
Rq = 8 kpc, and assuming p^ follows an NFW profile, 
when Ait!o = ±1 kpc, we find Awq ~ ±30 km/s and 
Apy/py Ki =f2 0%, similar to the range of v^ explored in 



Section IIIC Therefore for higher (lower) values of vq, 
the exclusion curves in Fig. Hj are shifted upwards (down- 
wards) by an additional ^ 20%. Comparing with Figs. [4] 
and [5] we see that this correction is of a similar size to 
what we found for Acr„/cr„ when independently varying 
Vo for the elastic and momentum dependent scattering 
case when M^ > 30 GeV, although the spread in the 
limits remains small. For light and inelastic dark matter, 
the change when varying vq independently is ~ 200%, 
therefore the spread in limits shown in Fig. |4] is overes- 
timated by ^ 10%, but still remains significant. 



IV. UNCERTAINTIES IN THE FORM OF THE 
DARK MATTER VELOCITY DISTRIBUTION 

In the previous section we investigated the effect 
of astrophysical uncertainties whilst always assuming a 
Maxwell-Boltzmann distribution for the dark matter ve- 
locity distribution. Dark matter numerical N-body sim- 
ulations have shown that a typical velocity distribution 
exhibits global (the overall shape) and local departures 
('bumpy' features) from a Maxwell-Boltzmann distribu- 
tion, and in this section we wish to investigate how these 
affect exclusion limits. 

We use data extracted from the Via Lactea II (VLII) 
[50] and GHALO [51] simulations available from [M]- To 



facilitate a more direct comparison between the VLII and 
GHALO simulations, we use the GHALOg data from ^39] 
in which the maximum circular speed in the GHALO sim- 
ulation has been rescaled to agree with that from VLII. 
As explained in f3S], the lack of baryons in the simula- 
tions means that the dispersion of the extracted velocity 
distribution is lower than expected from a galaxy which 
contains baryons and dark matter, because it is expected 
that the baryons will create a deeper gravitational poten- 
tial well, increasing the velocity dispersion. Since we are 
interested in departures from a Maxwellian distribution, 
we will compare the simulation limits with those from a 
Maxwell-Boltzmann distribution with the same disper- 
sion. This corresponds to using vq — 184 km/s in Eqn. 
[6] In the previous section we found that different as- 
trophysical parameters can lead to significant deviations 
in the exclusion limits, therefore we will match these as 
closely as possible so that we are just observing the ef- 
fects of the simulation distributions. Following [32] we 
take px = 0.3 GeV/cm^, vq = (10.0,220 + 5.25,7.17) 
km/s and v^sc = 550 km/s and 586 km/s for VLII and 
GHALOg respectively. 

With the simulation data, it is possible to examine 
the effects of local structure in the dark matter velocity 
distribution. This was done in |39| by extracting the ve- 
locity distribution for one hundred spheres of radii 1.5 
kpc and 1 kpc for VLII and GHALOg respectively, cen- 
tred at a distance of 8.5 kpc from the galactic centre. In 
Fig. [6] the thin bands show the exclusion curves from Via 
Lactea II (upper panels) and GHALOg (lower panels) as- 
suming elastic and momentum dependent scattering for 
aU one hundred spheres for CDMS Il-Si (yellow), CDMS 
II-Ge (blue), CRESST-II (green) and XENONIO (red). 
We have checked the robustness of these bands by com- 
paring the strongest and weakest limit with the sixteenth 
strongest and sixteenth weakest limit at each mass, and 
find that the limits are essentially the same. The dashed 
line shows the limits derived from a Maxwell-Boltzmann 
distribution with the astrophysical and dispersion param- 
eters above. We find that VLII and GHALOg produce 
similar limits, which are generally more stringent than 
those from a Maxwell-Boltzmann distribution. 

The inelastically scattering case has been studied in 
detail in [39] , therefore in Fig. [t] we only plot the limits 
for CRESST-II for 6 = 130 keV (for CDMS H-Ge and 
XENONIO we find a similar spread in the limits). To 
give an estimate of what the typical limits are, we have 
plotted the limits excluding five and sixteen maximum 
and minimum cross sections at each mass as the blue 
dashed and red dotted bands respectively. As we might 
expect, the inelastic case shows much more of a spread 
in the cross section compared to the elastic case, and the 
limits from VLII and GHALOg show some differences. 
Interesting features to note are the bumpy features in 
the GHALOg exclusion limits, and near threshold we find 
that some spheres can not set any limits for both VLII 
and GHALOg. The black dashed line shows the exclusion 
limits from a Maxwell-Boltzmann distribution which as 
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FIG. 6: (Colour online). The exclusion limits from the Via Lactea II (upper panels) and GHALOs (lower panels) simulations. 
The legend is the same as Fig. [4] The solid bands shows the range of limits from the velocity distribution of 100 random 
spheres centred at 8.5 kpc. For comparison, the dashed line show the limits from a Maxwell-Boltzmann velocity distribution 
with the same astrophysical parameters as found in the simulation halos. The Maxwell-Boltzmann distribution generally sets 
less constraining limits on (t„. 



in the elastic case, is typically more conservative than the 
limits from the simulations. 

We emphasise here that we have not been comparing 
the simulation limits with the SHM, in which the circu- 
lar speed enters Eqn. [6] but rather to a best fit Maxwell- 
Boltzmann distribution which has the same velocity dis- 
persion. Until simulations include baryonic effects, which 
will increase the dispersion, we urge caution in making 
direct comparisons to the SHM. The velocity distribu- 
tions from the simulations are interesting since they show 
that the overall shape of the limits can change, rather 
than causing a spread about a central value as we found 
when varying the astrophysical parameters (cf Figs, m 
and wh. Furthermore since different experiments vary 
by different amounts, when comparing exclusion curves 
at one experiment with preferred regions of parameter 
space at another (for example, comparing CDMS II-Ge 
to the DAMA/LIBRA or CoGeNT preferred region), one 
should bear in mind that different velocity distributions 
may affect the limits from different experiments by dif- 
fering amounts. 



CONCLUSIONS 



If we are to understand the nature of dark matter, it is 
vital that we have a full understanding of the astrophys- 
ical uncertainties affecting dark matter direct detection 
experiments. For spin-independent elastic and momen- 
tum dependent scattering with M^ > 50 GeV, we have 
shown that the exclusion limits are robust against vari- 
ations in the galactic escape velocity Vesc and the Sun's 
circular speed about the centre of the galaxy vq (right 
upper and left lower panels of Figs. [2] and |4]) and under 
realistic variations in the form of the velocity distribu- 
tion (middle and right panels of Fig. [6]), with uncertainty 
^10%. The major uncertainty in this mass range arises 
from the error in the local dark matter density (a factor 
of 2). 

In comparison, for lighter masses, we found that un- 
certainties in Vq and v^sc and the velocity distributions 
from the numerical simulations can shift the exclusion 
curves horizontally by ~ 1 GeV at masses M^ ~ 10 GeV 
(upper left panels of Figs. l2l l4]andl6|. Similarly for in- 
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FIG. 7: (Colour online). The exclusion limits for CRESST-II, assuming inelastic scattering from 100 random spheres from the 
Via Lactea II (left panel) and GHALOs (right panel) simulations. The dotted red and dashed blue lines indicate the range of 
limits when we exclude five and sixteen of the highest and lowest cross sections at each mass. For comparison, the long-dashed 
black line shows the limits from a Maxwell-Boltzmann distribution with the same astrophysical parameters and dispersion. At 
M^ ^ 60 GeV, we find that some of the spheres are not able to set any limits. As in Fig.|6J we find that the Maxwell-Boltzmann 
distribution generally sets less constraining limits on a„. 



elastically scattering dark matter, we found the vertical 
shift in the exclusion curves when varying vq or v^sc is 
large (^ 100%), with lighter target experiments such as 
CDMS II being particularly affected. Variations in the 
velocity distribution can also lead to significant changes, 
as we explicitly demonstrated in Fig. [7] for the CRESST- 
II experiment. 

Inelastic and light dark matter are particularly sensi- 
tive to astrophysical uncertainties because the minimum 
speed needed to scatter is just below the galactic escape 
speed, so experiments only sample the tail of the veloc- 
ity distribution. Even though the recoil spectrum for 
momentum dependent dark matter is different from elas- 
tically scattering dark matter (Fig. IT]) they both respond 
to astrophysical uncertainties in a similar fashion for ger- 
manium, xenon and tungsten targets, because Vmin re- 
mains far from Vesc- Therefore models which only sample 
the tail of the velocity distribution should carefully ex- 
amine the effect of astrophysical uncertainties on their 
limits. 
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Appendix A: Experimental Details 

In this Appendix we discuss the four experiments 
which we will use as our reference experiments, and 
highlight particular uncertainties associated with each. 
These experiments all publish unbinned data, so we use 
Yellin's 'pmax' method [S3] to calculate the limit on the 
dark matter-nucleon cross section, since it generally pro- 
vides a stronger constraint than unbinned methods for 
the case of an unknown background. It produces similar 
limits to the 'optimum interval' method used by CDMS 
II, the 'maximum gap' method used by CRESST-II and 
XENONIO (in their elastic dark matter analysis), and 
was later used by XENONIO (in their inelastic dark mat- 
ter reanalysis). 



1. CDMS II-Ge 

We initially consider the four published running peri- 
ods (runs 118, 119, 123-124 and 125-128) by the CDMS II 
collaboration for scattering off germanium in the energy 
range 10 — 100 keV [Tl [5M7T| . The first analysis saw one 
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event at 10.5 keV; the second, one event at 64 keV; the 
third, no events; and the fourth, notably saw two events 
at 12.3 keV and 15.5 keV. 

It was reahsed that there was an error in the pubhshed 
analysis of the first two runs (118 and 119) and a com- 
bined reanalysis was performed, which had redesigned 
cuts, and lowered the energy threshold from 10 keV to 
5 keV [75]. After the cuts were applied, it was found 
that the two events at 10.5 keV and 64 keV were no 
longer present, but two new events at 5.3 keV and 7.3 
keV passed all the cuts. These were not thought to 
be real dark matter events. The 5.3 keV event is very 
close to the energy threshold so was thought to be from 
background leakage, while the 7.3 keV event had an un- 
physical negative phonon delay, but was passed because 
of a pathology in the cut definitions. We follow CDMS 
by using this reanalysis when setting limits rather than 
the two published runs."* It should be noted that due 
to the lower energy threshold of these first two runs (5 
keV rather than 10 keV), CDMS II-Ge is able to place 
stronger limits on light dark matter than is often cal- 
culated. Similarly, the removal of the event at 64 keV 
leads to stronger limits being set on inelastic dark matter 
since this event lay near the peak of the recoil spectrum 
dR/dER [16]. 

In the analysis of the final run, an improved measure- 
ment of the germanium detector mass found a 9% de- 
crease in the effective exposure. In addition, a correction 
of the neutron systematic errors led to a 4.5% increase 
in the effective exposure [73]. We therefore multiply the 
published exposures (excluding the final run analysis) by 
0.955. 

The effective germanium exposure of the combined 
analysis of runs 118 and 119 is given in Figure 10.1 of 
[72]. We find a good fit to this curve for Er > 5 keV is 
given by: 

^^ = 80.1-0.31£;,,-225.52e-^«/2.09_74 4g -i=;„/20.o 
0.955 ' 

For the exposure of the 2008 analysis (runs 123-124), 
we use 0.995 x 397.8 kg-days and use the fit to the effi- 
ciency from [74] . 

For the 2009 analysis (runs 125-128), we use an expo- 
sure of 612 kg-days and find 

eff = 0.38-1.3x10-^ EB-0.17e-^''^'^-^^-2.96e-^''^^-^^ 

is a good fit to the efficiency.^ 



2. CDMS Il-Si 

The silicon data from the first two runs was also re- 
analysed and two events at 34.9 keV and 94.1 keV passed 



*We have checked the exclusion Umits without the unphysical event 

at 7.3 keV and the differences are small. 
^Figure can be found in the online supplementary material of [T]. 



all the cuts [72]. The 2008 analysis saw no events [75] . 
These runs had a low energy threshold of 5 keV and 7 
keV respectively, so again, are capable of setting strong 
limits at low masses. 

The effective silicon exposure of the combined analysis 
of runs 118 and 119 is given in Figure 10.1 of 72J. We 
find a good fit to this curve for 5 keV < Ejj < 20 keV is 
given by: 

Exp5.20 = 11-29 - 28.76-^"/^-^^ 

and for 20 keV < Er < 100 keV is given by: 

Exp2o-ioo = (0.55 - 0.63e~^«/i^-i)(0.061£;j^, + 29.5) 

For the exposure of the 2008 analysis, we use 53.47 
kg-days and fit to the efficiency from Fig. 5.13 of j75] . 



eff 





2Er/30 - 7/15 

0.0075£;i^ + 0.125 

0.09£:i^ - 1.03 

0.02£;i^ + 0.02 

0.4 



Er<7 keV 

7 keV <Er< 10 keV 

10 kcV < Er<U keV 

14 keV < Er < 15 keV 

15 keV < Er < 19 keV 
19 keV < Er < 100 keV 



3. CRESST-II 

CRESST-II have released data from two modules from 
the prototype phase [77] and two modules from the com- 
missioning phase [78j . Here we only include the seven 
events in the energy range 10 — 100 keV from the com- 
missioning phase, since there was a significant upgrade 
to the experiment after the prototype phase. Combining 
the data leads to significantly weaker limits [TO] . 



4. XENONIO 

Finally we consider the liquid xenon based experi- 
ment XENONIO, which combines scintillation and ionisa- 
tion measurements to infer the nucleus recoil energy Er. 
There have been two analyses performed; the first consid- 
ered only elastically scattering dark matter [75] and saw 
10 events in the energy range 4.5 — 26.9 keV; the second 
analysis, which used the same data but with redesigned 
cuts considered inelastically scattering dark matter [HO] 
and saw 13 events in an energy range extended to 75 keV. 
We use an energy resolution for XENONIO given in |74j . 

The quantity C^ff is needed to convert from the mea- 
sured scintillation energy (measured in keVee), to the 
nuclear recoil energy (measured in keV). Unfortunately 
Ceff is poorly determined and presents the largest sys- 
tematic uncertainty in the results of the XENONIO ex- 
periment. The original analysis, and the energies quoted 
above assume a constant value £eff ~ 0.19, however 
more recently, [8TH83] have determined that Ceff de- 
creases at low energies, but there is still some uncertainty 
in its value. Refs. [HIIIHI] found Ceff -- 0.15 for Er < 10 
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keV, while Ref. [53] found that it continues to decrease 
from Leff ^ 0.15 at En = 10 keV to Leff -^ 0.10 at 
E]^ ~ 5 keV. At M^ ~ 10 GeV, the exclusion limits from 
Ceff — 0.19 and Ceff from Manzur et al. differ by more 
than an order of magnitude [83; . 

When setting limits in this paper, we use the average 
value of Ceff from [5TH55] and use the elastic analysis to 
set exclusion limits for momentum dependent and elas- 
tically scattering dark matter, and the inelastic analysis 
when setting limits for inelastically scattering dark mat- 
ter. 

We end this section by mentioning another liquid 
xenon based experiment, ZEPLIN-III ,84,, which uses a 
value of Ceff which is smaller at low recoil energy than 
the three measurements mentioned above. Therefore in 
contrast to XENONIO, we would expect the ZEPLIN-III 
limits to be stronger at low masses if they were to use 
Ceff from Manzur et al. However we don't reproduce 
their limits here since they take into account background 
subtraction which is beyond our capabilities. 



by 



Appendix B: Analytic formula for C(^h) 



Here we present an analytic formula for Q{Eji)^ defined 




({Er)^ —f{v + ve) 



1 / 111 1 I l\ 

^ f e-" /"« ~ /3e-^"-/^" j V <Ve 



V > Ve 



(Bl) 



(B2) 



and /3 = or 1, depending on whether the exponential 
cutoff is desired. An analytic formula for C(^-r) with 
/? = was presented in [8 5) and we agree with their 
results in this limit. In this paper, we take /3 = 1. 

It is convenient to define Xesc = Vesc/vo, Xmin — 

Vmin/Vo and Xe = Ve/Vo, whcrC Vi = \vi\. 



The normalisation constant A^ is given by: 



A^ = ^3/2 3 
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